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We calculated the nuclear spin-lattice relaxation rate in the Mgi_a;ALB2 system and found 
that the orbital relaxation mechanism dominates over the dipolar and Fermi-contact mechanisms in 
MgB2 , whereas in AIB2 due to a smaller density of states and strong anisotropy of boron p orbitals 
the relaxation is completely determined by Fermi-contact interaction. The results for MgB2 are 
compared with existing experimental data, whereas our results for the doped alloy and AIB2 can 
be used in future experiments to verify the theoretical predictions about the electronic structure of 
this system. 
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The discovery of superconductivity in MgB20 has stimulated a significant interest in the electronic structure of 
• these intercalated boron compounds. For the superconductivity, one of the most important parameters is the 
value N(Ep) = N of the electronic density of states (DOS) at the Fermi level in the normal state. Several 
groups have perfor-nafid local density-functional (LDA) calculations of this quantity and the results scatter around 
0.7 states/(eV-f.u.)Ha. Experimentally N can be determined in many different ways but the nuclear spin-lattice relax- 
ed ation (NSLR) rate 1/Ti measurements represent an excellent opportunity to check experimentally not only the total N 
and its partial components but also their anisotropy, i.e. the distribution between the in-plane and out-of-plane p or- 
bitals. Experimental numbers have been reported in Ref. where values of TT\ = 180, 155 and 165 K-sec for MgB2 
on n B (fi = 2.689/ijv) were obtained. These authors interpreted the relaxatian-.rates in terms of dipolar and orbital 
' contributions due to the corresponding Korringa r ^io and the already knowncra theoretical total N. Below by using 
the general formulas for T\ in the hexagonal crystal] we evaluate the relaxation rates for MgB2 and AIB2 compounds 
using an LDA calculation. We will show that in MgB2 the orbital relaxation rate is about three times larger than 
i the dipolar and the Fermi-contact rates with all other contribution being much smaller, whereas the Fermi-contact 

mechanism is dominating in AIB2. Our calculated total rates are in fair agreement with those measuredoU, thus 
T — I \ indicating that the calculated N in MgB2 is basically correct, whereas in AIB2, due to an unusually high anisotropy 
^ • of p states, we predict a much smaller relaxation rate fully determined by Fermi-contact mechanism. 

I The standard technical details of the LMTO-ASA calculations are similar to those described in Refs. |^. In addition, 
we carefully checked the sensitivity of ASA results to the parameters of calculations. We used different exchange- 
correlation potentials and inputs with different radii of the B sphere (both with and without empty spheres). Our 
values of DOS for s and p orbitals Ni m for the B site in MgB2 and AIB2 are listed in Table I. One can see that in 
MgB2 all p orbitals on the B site have a sizeable DOS, while in AIB2 only p z orbital is significant with N px « 0.1N pz . 
The s component in AIB2 becomes relatively more important compared to MgB2 resulting in the dominance of the 
Fermi-contact NSLR in AIB2, as we will show below. In both materials the contribution of d states to NSLR is very 
small. As for the Mg site, the s component of N is the most important, and we expect that the NSLR on 25 Mg is 
controlled by the Fermi-contact mechanism. However, in this paper we will focus on the n B NSLR. 

To calculate T\ according to a general prescription one has to estimate the (r~ 3 )i expectation values for different 
orbitals and the electronic density at the nucleus if 2 . (0) / 47r. These parameters (properly normalized with the densities 
of states) are the largest source of uncertainty in our calculation due to the notable dependence of (r~ 3 ) p on the 
q \ chosen radius of the B sphere tb (from rs — 2.15ao to tb = 2Agq it decreases by ~11%). In our calculations we used 
the largest rs that were possible without a significant distortion of the band structure, 2.4ao for MgB2 and 2.1ao for 
AIB2. For these radii we have ((ao/r) 3 ) p = 1.11 in MgB2 and 1.37 in AIB2. Although due to the poor convergence 
these averages probably still have a considerable inaccuracy of about 10-15% (which is close to the experimental errox 
bar), we believe that they are the most accurate that may be obtained in ASA. For comparison, the atomic valueEl 



for ((ao/r) 3 ) p in B is 0.78. The electronic densities on the nucleus ag<^g(0)/47r for MgB2 and AIB2 were, respectively, 
2.68 and 3.02. 

The contributions to the n B relaxation rate for the polycrystalline sample calculated using the formulas of Ref. |?] 
are given in Table II. The in-plane and out-of-plane p orbitals in MgB2 have similar densities of states, and hence 
the relative magnitude of opbital and dipolar contributions to NSLR is similar to the 3/10 rule for p states in a cubic 
crystal described by ObataQ. The Fermi-contact contribution is also important and amounts to 30% of the orbital 
term. The contributions from the d partial waves to the dipole and orbital relaxation rates were small (at the order 
of 1%) due to the low diagonal and off-diagonal densities of states {Nd/N p ) 2 ~ 0.02 and (N p d/N p ) 2 ~ 0.05. The 
quadrupole contribution to NSLR is negligible due to a rather small 11 B quadrupole moment. 
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The values of T\ obtained in such manner correspond to the theoretical 'bare' N which does not include the exchange- 
correlation enhancement. We have estimated the effective Stoner exchange parameter to be I = AE/m — 1.7 eV 
from the splitting of the bands at the L point in the external magnetic field. The corresponding Stoner enhancement 
of the uniform susceptibility (which enters the Tf 1 ) can be writtenl as S = 3/ [(1 - NI) (3 - 21 N)] « (1 - IN)~ a , 
where a w 1.62 in the 3D case and in the 2D case a = 2. In our case due to the mixed 2D and 3D character of the 
bands it is not clear what value of a should be used in our simple estimation. We used a — 1.9, resulting in the 
enhancement of T-f 1 by approximately 60% (Table II). The total calculated relaxation rate 81T0 -4 K-sec for MgB2 
should be compared with the experimental rates of 56-64-10 (K-sec) -1 measured at temperatures slightly above 
TJaD. The fact that such simple estimate gives a faster relaxation compared to experiments, even without taking into 
account any other mechanisms of susceptibility enhancement, may serve as a manifestation of a possible importance 
of unusual effects resulting in the lowering of the effective N. 

We obtained a very different relation between the different mechanisms of NSLR in AIB2 where no traces of 
superconductivity have been found and where no NMR data are available. According to our theoretical estimation, 
due to the sharp decrease of the p component of N compared to MgB2, the orbital and dipolar contributions to 
NSLR become very small, and the total NSLR in AIB2 is completely dominated by the Fermi-contact mechanism. 
The resulting NSLR rate is more than two times smaller than in MgB2 (see Table II). Corresponding numbers for the 
anisotropy parameter at the Fermi level N px /N pz are 0.1 for AIB2 and 0.73 for MgB2. This specific feature of AIB2 
on this stage is completely theoretical prediction and more experimental information needed to check it. We also 
performed the rigid band calculations of T x in Mg 2 .Al 1 _ ;r B2. From Fig. 1 one can see how different mechanisms of 
NSLR are progressing as a function of doping. The sharp decrease of N in the 2D sheets of the Fermi surfacea leads 
to a corresponding lowering of all contributions to NSLR, and at the point of the structural transition we expect a 
very large T\ . We expect that the new NMR experiments for this alloy will be a crucial test of our understanding of 
the electronic structure of this system. Simultaneously with the nuclear quadrupole resonanse data on this material 
(which is related to the anisotropy of total charges on different p orbitals) the general picture of anisotropy of p 
orbitals in normal states can be build. 

The above calculations have been done for a polycrystalline material. Because single crystals are becoming available, 
we include our estimations of the anisotropy in the angular dependence of NSLR rateQ A + B sin 2 9. For MgB2 we 
obtained B/A « —0.06, so that the NSLR is nearly isotropic. 

In conclusion, we have performed LDA calculations for the NSLR rate in the Mg a: Ali_ :E B2 system. We find that 
the orbital mechanism of relaxation dominates over the spin-dipolar and Fermi-contact mechanisms in MgB 2 , because 
the boron p orbitals at the Fermi level are distributed nearly isotropically and have a large DOS. With the values 
of 0.7 st./eV for the bare DOS at the Fermi level for MgB2 and 1.7 eV for the effective Stoner exchange parameter, 
reasonable agreement is obtained with the experimental relaxation rates. We estimate the overall error of our T\ 
calculations in MgB2 as 15%. Strong anisotropy of p states at the Fermi level and different leading mechanism of 
relaxation in AIB2 compared to MgB2 is predicted. 

When this paper was completed we became aware that similar results for NSLR in MgB2 were independently 
obtained in Ref. [HJ where the T\ on 25 Mg and the Knight shifts were also computed, as well as the small core 
polarization term for n B. 

We are grateful to I. Mazin and F. Borsa for valuable discussions. This work was carried out at the Ames Laboratory, 
which is operated for the U.S. Department of Energy by Iowa State University under Contract No. W-7405-82. This 
work was supported by the Director for Energy Research, Office of Basic Energy Sciences of the U.S. Department of 
Energy. 



TABLE I. Partial DOS for s and p orbitals at B site, 10" 3 (eV-spin-atom)" 1 
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Contributions to (TTi)" 1 [10" 4 /(K sec)] 
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FIG. 1. Different contributions to the total B relaxation rate in Mgi_ a: ALjB2. Thick solid line, orbital; dashed, dipolar; 
thin solid, Fermi-contact. 
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